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ABSTRACT 

Interpretation of the Fe II(UV)/Mg II emission ratios from quasars has a 
major cosmological motivation. Both Fe and Mg are produced by short-lived 
massive stars. In addition, Fe is produced by accreting white dwarf supernovae 
somewhat after star formation begins. Therefore, we expect that the Fe/Mg 
ratio will gradually decrease with redshift. We have used data from the Sloan 
Digital Sky Survey to explore the dependence of the Fe II(UV)/Mg II ratio on 
redshift and on luminosity in the redshift range of 0.75 < z < 2.20, and we have 
used predictions from our 830-level model for the Fe II atom in photoionization 
calculations to interpret our findings. 

We have split the quasars into several groups based upon the value of their 
Fe II(UV)/Mg II emission ratios, and then checked to see how the fraction of 
quasars in each group varies with the increase of redshift. We next examined the 
luminosity dependence of the Fe II(UV)/Mg II ratio, and we found that beyond 
a threshold of Fe II(UV)/Mg II = 5, and M 250 o < -25 mag, the Fe II(UV)/Mg II 
ratio increases with luminosity, as predicted by our model. 

We interpret our observed variation of the Fe II(UV) /Mg II ratio with redshift 
as a result of the correlation of redshift with luminosity in a magnitude limited 
quasar sample. 
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1. Introduction 

Recently, the number of observational efforts focused on Fe II(UV) /Mg II emission ratio 
measurements has increased due to the potential use of these ratios to trace star formation 
history. Both Fe and Mg are produced by short lived massive stars (SN type II), and 
additional Fe is produced by accreting white dwarf supernovae (SN type la) sometime later. 
Thus, one expects that the Fe/Mg ratio is low at high redshift, and gradually increases with 
decreasing redshift, with the increase starting at a redshift of about 4, corresponding to an 
age for the Universe of 1.5 Gyr (Hamann and Ferland 1993; Yoshii et al. 1998), or even 
as early as a redshift of 6 if star formation in the extreme environment of a quasar central 
region begins at z > 10 (Matteucci & Recchi 2001). 

Although many independent ways to measure Fe II(UV) /Mg II ratios observed in the UV 
to IR range have been developed and applied to different data sets, the common conclusion 
is that the Fe II(UV)/Mg II ratio shows a large scatter at all redshifts, and little evolution 
with redshift. (Kwan & Krolik 1981; Wills et al. 1985; Kinney et al. 1991; Kawara et al. 
1996; Thompson et al. 1999; Iwamuro et al. 2002; Dietrich et al. 2002; Freudling et al. 2003; 
Barth et al. 2003; Dietrich et al. 2003; and Maiolino et al. 2003). Furthermore, it is usually 
assumed that there is a linear dependence between the Fe II(UV)/Mg II and Fe/Mg ratios 
because the ionization potentials of Mg (7.65 eV) and Fe (7.87 eV) are nearly the same. 

In our previous study, we showed that the Fe II(UV)/Mg II ratios are more sensitive 
to other physical properties of the emitting region than to abundance (Verner et al. 2003; 
2004). If the physical conditions are different in the Broad Line Regions (hereafter, BLR) of 
different quasars, the resulting scatter of Fe II(UV)/Mg II ratios obscures any dependence 
on abundance. Thus, it is important to explain the origin of the observed scatter before 
attempting to derive the Fe/Mg relative abundance in quasars. 

Although the prominent emission lines observed in quasars generally allow us to evaluate 
the metalicity of galactic nuclei, to link specifically the Fe II emission with Fe abundance 
is not a simple task. Compared to many other ions, the Fe II ion has a very rich spectrum 
due to its half-filled 3d-shell. As a result the Fe II(UV) band from 2200-3000A [hereafter, 
Fe II(UV))] can contain hundreds of strong lines. The Fe Il(UV)/Mg II emission ratio 
is therefore heavily affected by many Fe II lines but by only the two Mg II doublet lines 
at 2800A. Nevertheless, in order to measure the Fe/Mg abundance ratio, obtaining the 
Fe II(UV)/Mg II ratio is unavoidable. 

As quasar spectra (line intensities and continuum) are heavily affected by Fe II emission, 
it is possible to learn more about BLRs in quasars by studying how the Fe II originates. To 
achieve such a goal we have constructed an 830-level model for the Fe II atom and investigated 
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how Fe II(UV)/Mg II ratios vary with changes in hydrogen density, microturbulence and 
abundance. The model also predicts that the Fe II(UV)/Mg II ratio strongly depends upon 
the ionizing flux of the central source. 

In this Letter we have investigated if there is any dependence of Fe II(UV)/Mg II ratios 
with redshift. For the first time our approach combines model predictions with measurements 
[provided by Iwamuro (2004)] of the Fe II(UV)/Mg II ratios of quasars in the extended 
redshift range, 0.75 < z < 2.20, for quasars in the Sloan Digital Sky Survey (hereafter, 
SDSS) 1 . 

2. Fe II(UV) /Mg II Dependence on Redshift 

The Fe II spectrum has a rich and complicated mixture of forbidden and permitted lines. 
The pioneering work on the term analysis of singly ionized iron (Russell 1926) reported 61 
energy levels. Since work on Fe II energy measurements is still not completed (Johansson 
1978, 2004) it is not surprising that the very first attempts to explain Fe II in BLRs were 
not able to reproduce large Fe II(UV) /Mg II ratios when assuming solar abundance (Kwan 
& Krolik 1983; Wills, Netzer, & Wills 1985). 

Our 830-level model for Fe II atom in photoionization calculations is able to account 
for non- abundance factors (e.g. microturbulence, hydrogen density, and ionizing flux). This 
model is a natural extension (up to 14.1 eV) of our earlier model for Fe + , which included 371 
energy levels below 11.6 eV (Verner et al. 1999). The increase in the number of transitions 
from 68,638 to 344,035 is largely due to the increased density of energy levels at higher energy. 
The Fe II(UV)/Mg II ratios have been investigated assuming ranges of iron abundances (1, 
5, 10 in solar units) and microturbulence v tur b = 0, 1, 10, 100 km s -1 . Based on our model 
predictions we conclude that the Fe II(UV)/Mg II ratio is less sensitive to abundance than 
the Fe II(UV)/Fe Il(Opt) ratio (where Fe Il(Opt) is at 4000-6000A) (Verner et al. 2003, 
2004). 

After we found that the abundance is not the strongest factor that determines the 
Fe II(UV)/Mg II ratio in BLRs of quasars, calculations were repeated for a wide range of 
physical conditions, assuming solar abundance, and constant v tur b — 5 km s -1 , hydrogen 
density range n# = 10 9 5 — 10 13 '°cm~ 3 , total column density, Nh = 10 24 cm -2 , and the flux 
of hydrogen ionizing photons at the illuminated face 10 17 5 — 10 22 photons cm~ 2 s _1 . We 
employ the characteristic AGN continuum described in Korista et al. (1997), which consists 
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of a UV bump peaking near 44 eV, a /„ oc v~ l X-ray power law, and a UV to optical spectral 
index, a ox = —1.4. Figure 1 presents the results of model calculations of the Fe II(UV)/Mg II 
ratio as a function of density and ionizing photon flux. Although solar abundance is assumed 
throughout, our model predicts that the Fe II(UV)/Mg II may reach values as large as 40 
under conditions of high density and high ionizing photon flux. Most observed values lie 
in the range 2 < r < 5, with only a small fraction falling in the range 10 < r < 20. The 
Fe II(UV)/Mg II ratio should be relatively constant at low density and weak ionizing photon 
flux, but should rise steeply with ionizing photon flux beyond $ = 10 20 photons cm~ 2 s _1 
and n H > 10 n cm~ 3 . 

Our large Fe II model used in photoionization calculations predicts that the Fe II(UV) /Mg II 
ratios can have the same value over a wide range of physical conditions (Fig. 1). Therefore 
the typical observed Fe II(UV)/Mg II ratio of ~ 4 may arise from different excitation and 
density regimes. 

The SDSS database (Schneider et al. 2003) includes 11,677 objects within 0.75 < z < 
2.2. Iwamuro (2004) was able to measure the Fe Il(UV)/Mg II ratio for 10,670 of these 
quasars, and these are included in the sample under study. The smallest and the largest 
selected redshifts are z min = 0.7485 and z max = 2.1964, correspondingly. The median redshift 
of the sample is z me d = 1.4573. 

In Fig. 1, we see that at high flux levels, (log $ > 20.5), the Fe Il(UV)/Mg II ratio 
behaves differently in the high (log n H > 11) and low (log n H < 11) density regimes. At 
high density, collisions enhance the Fe II(UV) emission, while at low density, the Fe II(UV) 
emission drops relative to the Mg II emission. Very large values of the Fe II(UV)/Mg II 
ratio, those greater than 10, are expected only at high luminosity and high density. At lower 
flux levels (log $ < 20.5), the Fe II(UV)/Mg II ratio shows little density dependence, and 
we can generally conclude that the Fe II(UV)/Mg II ratio depends on the ionizing flux. 

We divide quasars into several groups depending on the value of ionizing flux, keeping 
the Fe II(UV)/Mg II ratios within a relatively small range: 1) Fe II(UV)/Mg II< 2 group 
shows almost no dependence on hydrogen density, and is mainly due to low ionizing flux; 2) 
2 < Fe II(UV)/Mg II < 5 group includes most quasars; 3) 5 < Fe II(UV)/Mg II < 7.5 group 
has a ratio a little bit higher than the typical value; 4) 7.5 < Fe II(UV) /Mg II < 10 ratios 
in this group in previous studies were usually interpreted as having increased Fe abundance; 
5) Fe II(UV)/Mg II > 10 group with large ratios that are relatively rare. 

The measured Fe II/Mg II ratio depends strongly on the fitting technique. High values 
of the Fe II/Mg II ratio result when the Fe II contribution to the Mg II doublet is taken 
into account using the quasar template from Wills, Netzer, Wills (1985). Results based 
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on the I ZW 1 template (Vestergaard & Wilkes 2001, Dietrich et al. 2002) underestimate 
contribution of Fe II to the Mg II doublet (see Verner et al. 2003; 2004 for more detailed 
discussion). 

Iwamuro et al. (2002), have developed a parameterized fitting procedure to obtain 
the Fe Il(UV)/Mg II ratio. The Fe Il(UV)/Mg II ratio is defined over 2150-3300 A to 
include the majority of Fe II(UV) emission and to avoid the feature below 2150 A. This 
feature is possibly due to the contribution of a Comptonized accretion disk (Zheng et al. 
1998). Iwamuro et al. applied their method to independently obtained data in the wide 
redshift range, < z < 5.3, [Kinney et al. 1991 and HST/FOS archival data (z < 0.17) 
, SDSS archival data (0.75 < z < 2.29), Thomson et al. 1999 data (3.1 < z < 4.7) and 
the high redshift sample (4.4 < z < 5.3)]. As a result of such examination they concluded 
that Fe Il(UV)/Mg II emission ratios in quasars show large scatter from 1 to 20 with little 
evidence of evolution of the Fe II(UV)/Mg II ratio with redshift, and that luminosity does 
not have a large effect on the Fe II(UV)/Mg II ratio, except for extremely bright objects. 

We have re-examined the redshift dependence of the Fe II(UV)/Mg II ratios measured 
by Iwamuro (2004) in the SDSS quasars over the range 0.75 < z < 2.20. We compute the 
fraction of quasars with various Fe II(UV)/Mg II ratios in each redshift bin (Az ~ 0.2). 
We find that the fractions are constant at low redshift, but show an increase in the relative 
numbers of quasars with Fe II(UV)/Mg II> 5 at redshifts greater than z ~ 1.8 (Fig. 2). 
This goes in the opposite direction of that predicted by the supernova enrichment scenario, 
which predicts that the excess Mg production in the early Universe would cause the Fe/Mg 
to decrease with redshift. 

In any survey which is magnitude limited, the most luminous objects will be seen over 
the entire survey volume, while the least luminous objects will only be seen nearby. This 
induces a correlation between luminosity and redshift. We adopt a cosmology with H = 
70 km s" 1 Mpc -1 and Q m = 0.3, Q\ = 0.7 for the estimation of the absolute AB magnitude at 
rest 2500 A, M 2 5oo- The smallest and the largest M2500 in the sample are M^qq = —21.67 mag 
and M^oo — —28.64 mag. The median luminosity of the sample is M^jo = —24.86 mag. 

For each Fe II(UV)/Mg II ratio, we plot the mean absolute luminosity of the quasars 
in each redshift bin vs. redshift, and find a very strong correlation between redshift and 
luminosity (Fig. 3). 

We next examine the luminosity dependence of the Fe II(UV)/Mg II ratio in the SDSS 
quasars over the range 0.75 < z < 2.20. We compute the fraction of quasars with various 
Fe II(UV)/Mg II ratios in each luminosity bin. We find the fractions are constant at low 
luminosity, but show an increase in the relative numbers of quasars with Fe II(UV)/Mg II > 
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5 at high luminosity, M 250 o < —25 mag. We have found that beyond a threshold, the 
Fe II(UV)/Mg II ratio increases with quasar luminosity (Fig. 4). This behavior is entirely 
consistent with our model (Fig. 1). At low flux levels, log$ < 19, increasing luminosity has 
little effect on the Fe II(UV) /Mg II ratio, while at higher flux levels, log$ > 19, the ratio 
steeply increases. 

The much larger number of atomic levels for Fe II atom used in the photoionization 
calculations provides an enhanced accuracy to account for the radiation field and enables 
us to explain the large Fe II(UV) /Mg II ratios better than in any previous efforts. The 
complicated level structure of Fe II ion leads to a more rapid increase of Fe II emission with 
increasing of ionizing flux. More Fe II levels become populated and more lines contribute to 
the total Fe II emission. As no iron overabundance is needed to explain the large Fe II(UV) 
/Mg II ratios that are observed, we have used solar abundances throughout our study (Verner 
et al. 2003, 2004). 

Observations of Fe II (Opt) are needed to constrain whether the discovered change of 
Fe II(UV) /Mg II ratio with luminosity also reflects a change in other physical conditions 
(e.g. hydrogen density) in emitting regions of BLRs. 

3. Conclusions 

While the extensive quantitative comparisons between Fe II(UV)/Mg II ratio measure- 
ments and model predictions are needed, it is clear that our model predictions are in a 
general agreement with our observational results. Although it has been suggested to use 
the Fe II(UV)/Mg II ratio to trace evolution, the physical processes forming lines in the 
Fe II(UV) band means that the Fe II(UV)/Mg II ratio tells us more about the central ion- 
izing source than about abundances. Thus, we can only convert an Fe II(UV)/Mg II ratio 
into an Fe/Mg ratio when a measurement of the Fe Il(Opt) band has also been obtained. 

Our analysis demonstrates that the observed Fe II spectra can be produced by a wide 
range of ionizing flux regimes. The statistical approach applied to physically distinguished 
Fe II(UV)/Mg II ratio groups will provide a more efficient way to improve our understanding 
of BLR in quasars than the universal template approach. 

Due to complexity and richness of the Fe II spectra in quasars, only the comparison 
between model predictions and extensive observations make it possible to provide accurate 
measurements of changes of physical conditions and abundances in quasars with redshift. It 
is possible that the change we see is due to not only luminosity but to the hydrogen density 
changes in the BLR with redshift as well. Wide wavelength coverage in the observations is 
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needed for a successful understanding. 

We find an increase in the Fe II(UV)/Mg II ratio with redshift, while evolutionary 
models predict a decrease. The Fe Il(UV)/Mg II ratio is not sensitive to abundance changes 
(see Fig. 4, Verner et al. 2003), but it is strongly affected by luminosity. We have found that 
the relative numbers of quasars with a high Fe II(UV)/Mg II ratio increases with luminosity, 
as predicted by our model. 

We conclude that the apparent change in the Fe II(UV)/Mg II ratio with redshift is 
simply a result of the correlation of redshift with luminosity in the magnitude limited quasar 
sample, and is not produced by a change in abundance. 

We are very grateful to Dr. Fumihide Iwamuro for providing us with his measurements 
of the Fe II/Mg II emission ratios in the spectra of SDSS quasars. 
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the central source affects the Fe II/Mg II emission ratio (Iwamuro et. al 2002). 
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Fig. I. — The Fe II(UV)/Mg II ratios predicted using our 830-level model for Fe + in a quasar 
BLR as a function of the hydrogen density % and flux of hydrogen-ionizing photons. The 
chemical abundances are solar, the cloud column density is 10 24 cm -3 , and v tur b = 5 km s _1 . 

Fig. 2. — The fraction of quasars in each redshift bin for various Fe II(UV)/Mg II ratio 
groups. The numbers on the top show the total number of quasars in each bin. The quasar 
sample comprises those quasars from from the SDSS in the redshift range 0.75 < z < 2.20. 

Fig. 3. — The mean absolute luminosity for each Fe II(UV)/Mg II group of the quasars 
in each redshift bin vs. redshift. The quasar sample comprises those quasars from from the 
SDSS in the redshift range 0.75 < z < 2.20. 

Fig. 4. — The fraction of quasars in each luminosity bin for various Fe II(UV) /Mg II ratio 
groups. The numbers on the top show the total number of quasars in each bin. The last 
bin with only 4 quasars has been ignored. The quasar sample comprises those quasars from 
from the SDSS in the redshift range 0.75 < z < 2.20. 
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